Giant fluctuations of superconducting order parameter in 
ferromagnet/superconductor single electron transistors 
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Spin dependent transport in a ferromagnet/superconductor/ferromagnet single electron transistor 
is studied theoretically with spin accumulation, spin relaxation, gap suppression, and charging effects 
taken into account. A strong dependence of the gap on the magnetic state of the outer electrodes is 
found, which gives rise to a negative magneto-resistance of up to 100 %. We predict that fluctuations 
of the spin accumulation due to tunneling of quasi-particles can play such an important role as to 
cause the island to fluctuate between the superconducting and normal state. 
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Electronic devices, which exploit the spin as well as 
charge of the electron have attracted much interest over 
the last several years. This interest in "spintronics" is 
stimulated by expected industrial applications, such as 
magnetic tunnel junctions based sensors, magnetic ran- 
dom access memory, as well as basic physics, offering new 
ways to manipulate quantum states using these two fun- 
damental properties of the electron. 

In nanoscale magnetic devices, such as single electron 
transistors (SET's), the presence of Coulomb blockade 
makes the charge transport sensitive to a very few elec- 
tron states, which results in a number of pronounced 
effects in magneto-conductance (see review [ij and refs. 
therein). The sensitivity to spin imbalance in the SET 
should be greatly enhanced if the island is made super- 
conducting, since the dependence of the gap on the num- 
ber of accumulated spins is highly non-linear. Takahashi 
et al. analysed a double tunnel barrier structure with 
ferromagnetic outer layers and superconducting central 
layer (F/S/F) without charging effects. They found that 
the anti-parallel alignment of the ferromagnetic films can 
produce a non-equilibrium spin distribution in the cen- 
tral layer, which suppresses the superconducting gap and 
thereby strongly influences the transport. Immamura et 
al. 01 considered the parity effect in a ferromagnetic SET 
with a superconducting island (F/S/F SET) in the limit- 
ing case of charging energy larger than the superconduct- 
ing gap, low voltage, and zero temperature. They found 
that the spin accumulation is zero in the limit considered, 
i.e., the superconducting gap of the island is unaffected. 
Magnetoresistance (MR) can exist due to the polarized 
spin of the single excess electron but is reduced in mag- 
nitude from the Julliere value by spin relaxation. 

In this letter we develop a comprehensive model of 
the F/S/F SET by taking into account spin accumula- 
tion, spin relaxation, gap suppression, charging effects as 
well as conduction fluctuations in experimentally realistic 
structures, at flnite temperatures. We flnd that spin ac- 



cumulation in the F/S/F SET suppresses the gap, which 
strongly affects the transport properties and results in 
a greatly enhanced MR. Progressively smaller spin accu- 
mulation is required to suppress the gap as the size of the 
island is reduced, making it feasible to electrically detect 
just a few non-equilibrium spins. We flnd that fluctua- 
tions of the spin accumulation can play such an important 
role as to cause giant fluctuations of the gap between the 
full value at a given temperature and zero (between the 
superconducting and normal state at fixed bias and gate 
voltages), with corresponding large fluctuations in the 
conductance of the device. Even in the parallel magnetic 
state of the device, where the average spin accumula- 
tion is zero, fluctuations of spin accumulation can lead 
to a signiflcant gap suppression. The transistor exhibits 
a practically full gate-controlled MR effect (MR~ 1, cor- 
responding to ~ 1000% in typical GMR notations). We 
compare the results of the model with the recent, rather 
divergent experimental flndings on ferromagnetic SET's 

HllQ. 

The model geometry is shown in Fig. 1(a). Two ferro- 
magnetic electrodes are connected to a superconducting 
island through tunnel junctions. The tunnel resistance 
is assumed to be much larger than the Quantum resis- 
tance, i?, > i?Q = h/e^ ~ 26 kn. We assume collinear 
magnetization of the two electrodes, the magnetization 
of the left electrode can be either parallel (P) or anti- 
parallel (AP) to the magnetization of the right electrode. 
The resistance of the tunnel junctions is spin dependent: 
R] = 2RJ{1 + P,), r\ = 2i?,/(l-P,) (where i = L (left) 
or R (right) denotes the junction). Here the polarization 
is of the same sign {PlPr > 0) for the parallel ahgn- 
ment. Letting Pl —Pl (or Pr —Pr) for reversed 
magnetization of the left (or right) electrode, we obtain 
the resistances in the anti-parallel alignment. 

The spin relaxation time can be comparable to or 
greater than the average time between tunnel events, so 
that a non-equilibrium spin distribution can accumulate 
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on the island, resulting in a shift in the chemical po- 
tential. The shift in chemical potential is schematically 
shown in Fig. 1(b). We consider the regime where the 
energy relaxation time is much shorter than the average 
time between two successive tunnel events (no gap sup- 
pression is expected in the opposite, elastic regime; for 
details see 0) • The electrons are then instantaneously 
thermaHzed after they have tunneled and are described 
by the Fermi distribution function. This shift in chemical 
potential, A/i'^, is related to the deviation in the number 
of spins from the equilibrium value, n°', for each spin- 
band, a |: 

/oo 
dEN{E)[f{E-Afi'')-f{E)], (1) 
-OO 

where N{E) = e{E - A) ^ Jfj^, is the BCS density of 
states, p is the density of states at the Fermi energy (per 
spin) in the normal state, f{E) is the Fermi distribution 
function, V is the island volume. The shift in the chemi- 
cal potential on the superconducting island influences the 
gap through the BCS gap equation 

In (^) = ^ [m + Am^) + m + , 

(2) 



where ^ = y/E^ + A'^ and is the Debye frequency. 
For arbitrary spin distribution on the island, (arbitrary 
andn^), solving equations (QJ and l|2l self-consistently 
yields the shift in the chemical potential, Afj.'^ , and the 
gap, A. In Fig. 2 the chemical potential and the gap as a 
function of the spin imbalance, s = — n^, are shown for 
the case rt^ — —nK When the spin imbalance increases, 
there will be a critical value of s where the gap will dis- 
continuously drop to zero and the chemical potential will 
drop to a lower value. After the complete suppression 
of the gap the chemical potential will increase linearly 
with the slope of l/pV as in the normal conducting case. 
Reducing the spin imbalance will restore the gap at the 
same critical value of s. Note that the gap and chemi- 
cal potential do not only depend on and through 
the combination s ~ — , the total spin accumula- 
tion. Different and having the same s can produce 
different A and Ap'^ . 

The transport properties are determined bvthe tunnel- 
ing rates according to the "orthodox theory" [Sj . The tun- 
nehng rates depend on the charge of the island through 
n = +n-^, as well as and separately, through the 
dependence of A and Ap"' on and nK The tunneling 
rates are spin dependent also through the spin dependent 
resistances: 
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where 5E°^^°^{n) is the change in the charging energy 
for tunneling one electron on or off the island through 
the left or right junction: 
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Here Ce ^ Cl + Cb. + Cg is the total capacitance of the 
island and Ug — CgVg/e is the normalized gate charge. 
Tunneling that increases the the energy have exponen- 
tially small rates (zero at T = 0) , giving the characteris- 
tic Coulomb blockade below a threshold voltage. 

In what follows we choose the parameter space to rep- 
resent the typical nano-lithographically fabricated SET 
structures. We therefore consider SET's with island di- 
mensions in the 10^ — 10^ nm'^ range, such that up to a 
hundred states will participate in the transport. The rel- 
atively large number of states involved makes the Monte 
Carlo method ^3 more suitable than the master equa- 
tion approach. The intertwined dependence of the gap. 



the chemical potential and the spin accumulation can be 
included in a straight forward way with the Monte Carlo 
method. Moreover, the fundamental stochastic nature 
of tunneling and the resultant fluctuations in the spin 
transport through the structure are well modeled by the 
Monte Carlo procedure. 

The state of the system is determined by -n) and n^, 
from which the gap and shift in chemical potential are 
calculated using Eqs. and After obtaining A and 
Ap'' , the eight possible tunneling rates are computed, 
which increase or decrease or by one. The inter- 
val zero to one is then divided into eight segments with 
widths proportional to the tunneling rates and a random 
number in the interval [0,1] is generated to choose one of 
the eight tunnel events. The probability to preserve the 
present configuration decays as P = exp[— (i — io)rs], 
where to is the time of the previous tunnel event and Fs 
is the sum of the eight tunneling rates. A random num- 
ber in the interval [0,1] gives the time between successive 
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tunnel events by inverting the equation for the proba- 
bility. The spin flip process is then accounted for in a 
phenomenological way by letting s sexp[— (t~to)/'^sf], 
where Tgf is the spin flip time. This procedure is repeated, 
keeping track of all quantities, until good statistics is ob- 
tained. 

The dynamic effect of gap suppression is shown in 
Fig. 3, where the spin accumulation and superconducting 
gap are plotted versus time at flxed bias and gate volt- 
ages. As soon as the spin imbalance reaches the critical 
value (~ 18 for these device parameters) the gap drops 
to zero and the island goes normal. At the critical bias, 
fluctuations of only a few polarized spins on the island 
cause giant fluctuations of the order parameter from the 
full value at the given temperature to zero. This effect 
should be directly observable in experiments employing 
time resolved (<l/is scale) conductance measurements. 
Typical current- voltage (I-V) measurements (>1 ms time 
constant) will record the gap, which is a time average of 
the dynamic gap shown in Fig. 3. Fig. 4 shows this mean 
value of the superconducting gap, (A), as well as its stan- 
dard deviation, <5A = [(A^) - (A)^]^/^, as a function of 
bias voltage. The fluctuations result in a smooth rather 
than discontinuous transition from the superconducting 
to the normal state as the spin accumulation is increased 
in the AP state by increasing the bias. The fluctuations 
peak at the critical bias in the AP state (0.42 mV in this 
case). Contrary to what would be intuitively expected, 
fluctuations of the spin accumulation on the island in 
the P state can be very large. The fluctuations increase 
with bias voltage and becomes as large as the critical spin 
accumulation giving a signiflcant gap suppression. This 
unexpected effect is a consequence of the stochastic dy- 
namics of spm-dependent transport through the device 
(charge fluctuation is suppressed by the charging energy) 
and the highly non-linear dependence of the gap on spin 
accumulation, since the time averaged spin accumulation 
in the parallel state is zero. 

The suppression of the superconducting gap by spin 
imbalance in the AP state results in fluctuations of the 
conductance of the device, which are greater in magni- 
tude as the conductance itself. Fig. 5a shows the time 
averaged I-V characteristics of the ferromagnetic SET 
in both the AP and P states. At the critical bias, the 
gap suppression leads to a 3-fold increase in current as 
a function of the magnetic state of the device, which 
results in a correspondingly large negative magnetore- 
sistance. An even stronger MR effect, MR>0.8 is ob- 
tained at the critical bias by applying a gate voltage as 
shown in Fig. 5b. The island size chosen for the illustra- 
tion above is representative of typical experimental im- 
plementations today, but not the smallest fabricated to 
date. For smaller superconducting islands, progressively 
fewer spins are needed to produce the full effect of gap 
suppression, with MR reaching « 1. Increasing the spin 
relaxation time or reducing thermal excitations by going 



to lower temperatures enhances the effect. 

To estimate the parameters of a realistic experimen- 
tal realisation of a F/S/F SET, assume a SET fabricated 
with Cobalt leads having a spin polarization of 0.35, Alu- 
minum as the superconducting island with A = 200 /ieV, 
p = 10^^ [m^eV]~^ and a volume of 15 x 70 x 200 nm^, so 
that pV = 2.1 • 10^ [eV]"^ The spin relaxation time of 
quasi-particles in superconducting Al has not been mea- 
sured, but for normal Aluminum at 4.2 K spin relaxation 
times between 10 ns and recently 0.1 ns have 
been reported. Assuming a spin relaxation time of Tgf « 
1 ns, T =0.05 K, i? = 30 kf], and = 100 ^eV this 
sample would have the transport characteristics shown 
in Figs. 2-5. If the spin relaxation time was shorter, the 
island volume would have to be smaller to achieve the 
same effect. Increased Tgf in the sup erconducting state 
has been predicted theoretically |lj|, which would lead 
to stronger gap suppression and higher MR. 

When we compare our theoretical results with the re- 
cent experiments on F/S/F SET's we draw the con- 
clusion that the reported MR is unlikely due to gap sup- 
pression by spin accumulation. The volume of the islands 
studied was rather large in both experiments {pV « 10^ 
[eV]~^) and the temperature relatively high (T « 0.25 
K). The combined effect of large volume and high tem- 
perature would give insigniflcant gap suppression in the 
voltage interval where large MR was reported. There are 
several other points where our simulations do not agree 
with the interpretation of experiments reported. In both 
experiments the maximum MR was observed at or below 
the threshold voltage, eVth < 2 A -t- -Ec- In our calcula- 
tions we always flnd Vth to be the same in the two states, 
and the maximum MR occurs above Vth- This arises 
form the fact that the spin accumulation is directly pro- 
portional to the current, so there is no mechanism for 
spin accumulation below the threshold voltage, where no 
current flows. Further, the extracted gap in jj| did not 
go to zero monotonically with increasing voltage but had 
a maximum at an intermediate voltage, contrary to what 
we predict. The temperature dependence of the MR mea- 
sured in 0, with the peak of the MR occurring at the 
same bias voltage, is in contrast with our calculations 
where the MR peak shifts to higher bias with increasing 
temperature (not shown). More experiments are needed 
to resolve these controversies. 
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Figure 2: The self-consistent solutions of Eq. Q and 
The gap. A, and chemical potential, Afi, as a function of the 
spin imbalance for pV = 2 ■ 10^ [eV]"^ and T = 0.05 K. The 
critical spin accumulation, 5*0, is shown. 
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Figure 1: a) The single electron transistor consisting of two 
ferromagnet (F) electrodes and a superconducting (S) island, 
b) Schematic density of states for the ferromagnet (left and 
right) and superconductor (middle) for the anti-parallel align- 
ment of the electrodes, when Vl = V/2 and Vr = —V/2. 
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Figure 3: Time trace of the gap (solid line) and spin accu- 
mulation (dots) at in Fig. 2, corresponding to V =0.42 mV 
and Vg =0 mV. The device parameters are Rl = Rr = 30 
kn, Cl = Cr = 0.8 fF, Pl = Pr = 0.35, pV = 2 ■ 10® 
[eV]"^ T = 0.05 K, Tsf = 1 ns. Sc is indicated with a dashed 
horizontal line. 



0.2 



0.15 



0.05 



















— <AV 








— <A),p 








— - SAp 








--■5\P 


1 
1 
1 
I 








1 / 

1 ,' 






1 
1 
1 
1 
1 
I 


\ 1 

V 








1,1 



0.2 0.4 0.6 0.8 1 

V(mV) 



Figure 4: The time averaged value (solid) and standard devi- 
ation (dashed) of the superconducting gap for the P and AP 
states of the device in Fig. 3, VgCg/e = 0.5. 
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Figure 5: a) I-V characteristics in the P and AP state and 
MR of F/S/F SET with same parameters as in Fig. 3 and 
= 0; b) the gate dependence at V=0.45 mV. 



